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Abstract 


t 


Hnpiti developments in the fabrication technology have resulted in the 
scaling down of MOSFET device dimensions to ultra-small feature sizes. Thus, 
MOSFET models which accurately incorporate the various short-channel ef- 
fects, are becoming important. The abnormal increase observed in the transcon- 
ductance of ultra-short channel length MOSFETs is one such phenomenon 
which needs to be modeled accurately, as the transconductance parameter of 
the MOSFKT is extremely important from analog as well as digital circuit 
point of view. Previous authors have modeled this abnormal increase in the 
tmnsconduetance by assuming that there is an overshoot of the electron ve- 
locity beyond the saturation value in Silicon for ultra-short channel length 
devices. However, in their models, they assumed that the mobility is pinned 
to its low-field value, which is highly inaccurate, since the electric field present 
in the channel for ultra-short channel length MOSFETs is very large, and the 
assumption of the low-field mobility to be applicable for such high channel 
fields may not hold true. In this work, the abnormal increase in the transcon- 
ductance with decreasing channel lengths has been accounted for by proposing 
a new expression for the electron mobility, which is a function of the critical 
electric field, the channel length, the oxide thickness, the applied drain-to- 
source voltage, and the applied gate-to-source voltage. The results obtained 



from our model have boon compared with those reported experimentally, and 
a i\uud match between the two is seen. 

Another effect, which becomes significant, with shrinking device dimensions 
is the substrate leakage current caused by the impact ionization of high energy 
carriers in the surface depletion region near the drain. Previous authors have 
modeled the ionization length near the drain in a totally empirical manner. In 
this work, a simple and accurate approximation of the ionization length, based 
on a calculation of the electric field distribution near the drain region has been 
developed. Our model results are compared with the experimental results 
reported in literature and a good correlation between the two is obtained. 
This improved model for the ionization length has been used in the substrate 
current model developed in this work. A comparison of our simulation results 
with the ones reported experimentally shows a good match between the two. 
Prior to modeling these ultra-short channel effects, an existing physics based 
MOSFKT model available in literature is simulated in order to gain a better 
understanding of MOSFET device physics and simulation problems. 
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Chapter 1 


Introduction 


Wry Large Scale Integration (VLSI) is the assembly of millions of tran- 
sistors on a single chip known as the Integrated Circuit (IC). One of the main 
reasons for the popularity of the VLSI systems is the advancements in the fabri- 
cation technology which resulted in a very high packing density, made possible 
by a reduction in the sizes of the individual transistors on the VLSI chip. The 
metal-oxide-semiconductor field-effect transistor (MOSFET) has emerged as 
the most important electronic device, superseding the bipolar junction tran- 
sistor in sales volume and applications. Its strength derives from its simple 
structure, low fabrication cost, and low power dissipation. For these reasons, 
MOSFETs continue to be the most popular device for VLSI. 

Computer Aided Design (CAD) and simulation of electronic circuits are the 
main factors contributing to the success of VLSI. For accurate simulation of the 
electronic circuits, the devices used in the circuits must be modeled accurately. 
Thus, device modeling becomes an essential step in computer aided design and 
simulation of electronic circuits. 



1 h<* device model is a mathematical representation of the physical charac- 
teristics of the device. It is a set of mathematical equations which describe 
the physical behavior of the device under various input conditions. The set 
of parameters (constants) in these mathematical equations that represent the 
device phenomenon analytically are called the physical parameters, whereas 
the parameters which are used in the model equations such that these equa- 
tions are the numerical approximations of the device behavior are known as 
the empirical parameters. 

Development of models, which are totally based on the physical phe- 
nomenon, is a difficult task and also computationally very intensive. Although 
such models are highly accurate, they prove to be inefficient for circuit simu- 
lation. The totally empirical models, on the other hand, are computationally 
very efficient, however, they are less accurate and do not always give a true 
physical picture of events taking place within the device. With VLSI gaining 
popularity ami device dimensions continuously shrinking, the existing MOS- 
FET models, which hold their validity for long channel lengths, become grossly 
inadequate in properly simulating the device operation for submicron channel 
lengths. Therefore, some modifications of the existing models are necessary in 
order to accurately portray the electrical characteristics for such ultra-short 
devices. Some of the second order effects which are dominant in short channel 
devices are described below. 

In short-channel devices, it is found that in the saturation region, there is 
an appreciable increase in the drain current Ip with an increase in the drain- 
to-souree voltage Vos, particularly so for very short channel length devices. 
This phenomenon is termed as the channel length modulation effect, and is 
caused by a variation in the effective channel length due to the applied \ 'os. 
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'I his effort is either negligible or absent in long channel devices. 

In the subthreshold region, the lowering of the source-channel potential bar- 
rier by the applied drain-to-source voltage in short channel devices is known 
as the Drain Induced Barrier Lowering (DIBL) effect. The immediate man- 
ifestation of this effect is an increase in the subthreshold current. In effect, 
the threshold voltage Vp becomes a function of the channel length L and the 
applied with IV reducing with a reduction in L and with an increase in 
Yds- This phenomenon is totally absent in long channel devices. 

Another important effect that comes into prominence with decreasing chan- 
nel length is an abnormal increase in the device transconductance g m for 
ultra-small channel length devices [1-4]. The reason given for this higher than 
theoretically predicted maximum values of g m for ultra-small channel length 
devices, is the electron velocity overshoot [1-4]. In the model developed by 
Roldan ct a l. (5], an electron velocity overshoot parameter A a is used in order 
to model the carrier drift velocity v in the inhomogeneous electric field present 
in the channel. In the expression for the transconductance, A a accounts for the 
increase in the transconductance for ultra-small channel lengths, while assum- 
ing a constant value for the low-field electron mobility. This approximation 
of using the values for the low-field electron mobility for such high values of 
electric field present in the channel is grossly inaccurate. 

Another effect which becomes significant when device dimensions are scaled 
down is the substrate leakage current caused by the impact ionization by high 
energy carriers. When a MOSFET is biased in the saturation region, a mea- 
surable substrate current J 4U & always flows. This current is the result of the 
electron-hole pair generation in the surface depletion region near the drain due 
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to impact ionization. It is extremely sensitive to device dimensions, channel 
doping, terminal voltages, and temperature. In the substrate current model 
proposed by Arora and Sharma [6], the carrier ionization length U near the 
drain is modeled in a semi-empirical manner. Since the substrate current I su b 
depends exponentially on this ionization length Id, a small error in its calcula- 
tion can lead to a large error in the computed value of J 3U {,. As Id is modeled 
semi-empiricaliy, there is no physical significance for the model parameters, Iq, 

1 1 , and h appearing in (6). 

In this work, we have made an attempt to model the last two effects men- 
tioned in the previous paragraphs, i.e., the transconductance overshoot and 
the substrate current in ultra-short channel length devices, in an analytical 
manner, based on device physics. However, since totally physical models are 
difficult to develop and are computationally inefficient, fitting parameters are 
used whenever required. A brief description of the work is given below. 

The increase in the transconductance g m with decreasing channel lengths 
has been accounted for in this work by proposing a new expression for the 
electron mobility, which is a function of the critical electric field, the chan- 
nel length, the oxide thickness, the applied drain- to-source voltage Vos, and 
the applied gate-to-source voltage Vfej. With large channel electric field and 
small channel length, the mobility is bound to increase, and the increase in 
the transconductance is then accounted for by the correspondingly enhanced 
electron mobility. Our model results are compared with those reported in 
literature, and they show a good match. 

For the other effect, i.e., the substrate current due to impact ionization of 
high energy carriers, a simple and accurate approximation of the ionization 
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h'ngtii /,*, bast'd on a calculation of the electric field distribution near the 
drain region has been developed in this work. Our model results are compared 
with the MINIMUS simulation results reported by Wong and Poon [7], and 
a good correlation between the two is obtained. This new model for Id is 
then used to calculate the substrate current. The results obtained from the 
improved substrate current model are then compared with the experimental 
data reported by Arora and Sharma [6], and a good match is obtained between 
tin 1 two. 

However, before attempting to model these two effects for ultra-small chan- 
nel length MQSFETs, an existing physics based MOSFET model, reported by 
Joardar et al. [8] is rederived and simulated in this work in order to gain a bet- 
ter feel for the MOSFET device physics and the simulation problems. A major 
drawback in the existing MOSFET models is the presence of discontinuities in 
the derivatives of the drain current with respect to the gate-to-source voltage 
anil the drain-to-source voltage. These discontinuities arise primarily due to 
the nature of the velocity-field relations and the smoothing functions used in 
these models for the transition region between the linear and the saturation 
regions. The model developed by Joardar et al. [8] attempts to overcome 
t hese problems. The new model [8] is normalized surface potential based and 
is quasi-static in nature, i.e., the variation of the terminal voltages is assumed 
to be sufficiently slow, so that the they can be assumed to be identical to the 
tic voltages at any given time. This model is simulated in our work and the 
results obtained are compared with the experimental values given by Joardar 
et al. [8j. 

The organization of this report is as follows. The rederivation and simula- 
tion of the physics based MOSFET model developed by Joardar et al. [8], and 
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tin* comparison of the* simulated results with those reported [8] are presented 
in Chapter 2. In Chapter 3, the improved model for the ionization length la 
is given, which is used to obtain the improved substrate current model. This 
chapter also includes the comparison of the simulated model results with the 
experimental ones reported by Arora and Sharma [Gj. The enhanced mobility 
model developed in this work is presented in Chapter 4. This chapter also 
includes the comparison of the transconductance values obtained from the ex- 
pression derived in this work with the experimental transconductance values 
reported by Sai-Halasz et al. [1], Finally, Chapter 5 includes the summary of 
the entire work. The shortcomings of the models developed in this work along 
with the future scope for improvement, are also discussed in this chapter. 
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Chapter 2 


A Physics Based MOSFET 
Model 


2.1 Introduction 

There have been many MOSFET models developed over the years. Most of 
these models have drawbacks stemming from simplifying assumptions made 
during their development. These drawbacks limit the accuracy of these models 
in specific regions of device operation. A major drawback is the presence 
of discontinuities in the derivatives of the drain current with respect to the 
gate-to-source voltage and the drain-to-source voltage of the MOSFET. These 
discontinuities arise primarily due to the nature of the velocity-field expressions 
and the smoothing functions used in the transition area between the linear 
and the saturation regions in these models. A compact model is presented by 
Joardar H al {8], which overcomes the problems mentioned above. The new 



new model is based on the normalized surface potential and is quasi-static in 
nature, i.e., the variation of the terminal voltages is assumed to be sufficiently 
slow, so that they can be assumed to be identical to the dc voltages at any 
given time. 

In this chapter, the model developed by Joardar et al. [8] is rederived 
in order to get a better understanding of MOSFET device physics and the 
problems associated with device simulations. The simulated results obtained 
from I Ik* model are then compared with the experimental values reported by 
Joardar et al. [8j. 

2.2 The Static Drain Current Model 

Consider an n-channel MOSFET with the structure as shown in Fig.2.1. The 
channel is divided into three sub-regions, namely A, B, and C, with lengths 
l $ atsi Itmi and haw respectively. Let V D , Vo, Vs, and V B be the drain voltage, 
gate voltage, source voltage, and body voltage respectively. The regions A 
and C do not exist simultaneously, i.e., when Vjp > Vs, hats is zero and region 
A is absent; and when Vs > V D , l satD is zero and region C is absent; and 
for V D = Vs, both regions A and C vanish completely. The drain current I D 
saturates as Vp increases above a certain threshold value V Dsat . Therefore, in 
either of the regions A or C, the carriers are assumed to flow with the saturation 
velocity Vgat, depending upon whether Vs > Vo or Vs < Vd respectively. For 
Vd < V Djat , Id increases almost linearly with increasing V D (for small values 
of V 0 ). Thus, in region B, the carrier velocity is assumed to be less than v 3at . 

The general expression for the drain current Ip can be written in terms of 
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Figure 2.1: The details of a MOSFET structure showing the three regions A, 
B, ami 0, with / S0 ,s, h tn , and being their respective lengths. 

the drift and the diffusion components respectively (8] 


I D = -lFuQ'+ IV D 


d _9± 

dy ’ 


(2.1) 


where H' is the effective channel width, v is the average carrier drift velocity, 
Q[ represents the inversion charge density per unit area in the channel, and D 
is an effective diffusion constant. The normalized potential is defined as [8] 


Fi(y) = Q;(y)/c; 


ox> 


( 2 . 2 ) 


where C' t is the oxide capacitance per unit area and is given by 


C„ = <«/ t ox , (2.3) 

with a»d t ox being the oxide permittivity and thickness respectively. Thus, 
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the current equation can be rewritten as [8] 


b = -U'CnFr + WC„D^p-. (2.4) 

dy 

The expression for the carrier velocity as a function of the lateral electric field 
can be given as (8j 


Hejf(d4/dy) 

o + <S“g) 2 )‘' 2 ’ 


(2.5) 


where 4>{y) is the surface potential and d<p/dy is the lateral electric field. The 
effective mobility fi t fj can be given by [8] 


AW/ = £V/„, (2.6) 

where Uq is the low field mobility and / M is the mobility degradation factor due 
to the vertical (gate) field. Equation (2.5) satisfies the necessary conditions 
for a sufficiently accurate velocity-field expression, and also represents the 
measured characteristics adequately [9]. The mobility degradation factor fp is 
typically expressed in terms of the excess gate voltage beyond the threshold 
voltage. However, since this model [8] is surface potential based, fp can be 
expressed in a more fundamental form as a function of the vertical electric 
field, and is given by the empirical relation [8] 

f M = 1 + UBRED(E av x 10 -8 ) fiAVEA/> , (2.7) 

where U Bit ED and EAVEXP are dimensionless model parameters, and E av 
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is t lit* effective vertical electric field given by the expression [8] 


E«v = { Q' b + EAVFC^esi, (2.8) 

where ICAVF is another dimensionless fitting parameter, Q' b is the depletion 
charge per unit area, and is the permittivity of Si. The form of this model 
for the vertical electric field induced mobility degradation is suggested by Sun 
and I Mummer [10]. 

Following well known principles of MOS physics, the normalized potential 
F, can he expressed as [8] 


W) = I'cb - V FB - <f> - Q' b (<f>)/C' ox , (2.9) 

where O is the surface potential, Vqb is the applied gate-to-bulk voltage, and 
Yyn is the Hatband voltage. In the linear region, i.e., V Dsai > V D > Kg, the 
regions ,4 and C in Fig.2.1 are absent, and l sa ts = hato = 0. At y = 0 (i.e., at 
source), the exact values for the surface potential and the depletion charge per 
unit area are denoted by <f>, and Q' bS respectively. In the region 0 < y < l tin, 
Q b can be approximately given by [8] 

Q'b & (Qbs + Q&b)/ 2 + ^(^ — ^s/2 - ^d/2), (2.10) 

when* Q' b]) ami <pp are the depletion charge per unit area and the surface 
potential respectively at the drain end (i.e., at y = h in ), g is the dimensionless 
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(lUinmH's liiiouiizntion factor and a (C/V-cm 2 ) is given as [11] 


2 ‘ dip + <|0 /■ ( 2 - u ) 

Equal ion(2. 10) is derived in the following manner. First, Q' b (<t>) is linearized 
by expanding it in Taylor series around <f> = <f> s , and retaining only the first 
two terms. It is then similarly linearized around <f> = <j>o. Finally, these two 
linearized versions of Q' b are averaged in order to yield the given expression 
[ivqn.(2. 10)]. Thus, the expression for the normalized potential F* linearized 
in terms of the surface potential 4> can be given as [8] 


F<{4>) « (Fis + Fj £>)/2 - (1 + ga)((j> - <£ s /2 - ^ c / 2 ), ( 2 . 12 ) 

where F,s and F t p are the exact values of Fi at the source and the drain 
respectively. Substituting these values in the expression for the drain current 
Jo (Eqn.(2.1)] and integrating, one obtains [8] 


, |(F« + F iD )(p D - fe)/2 + V, h (F<s - fio)] 

where V t h is the thermal voltage, given by V th = kT/q, with k being the 
Boltzmann’s constant, T being the absolute temperature, and q being the 
electronic charge, L rt d is the channel length reduction factor, given by L re d = 
with L being the total channel length, and l lin bounding the ohmic 
region of the channel, as shown in Fig.2.1. In Eqn.(2.13), the effective mobility 
//,// is replaced by another term n f /f avg , which is an average of its values at 
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?/ - 0 and at y = /(,„ (given by ft r j}s and {i r /jD respectively), and is given by 

18 ] 


t l effavg — [t‘ef/S + Vc/Jd)/ 2- (2-14) 

The first and the second terms in the numerator of Eqn.(2.13) are the drift 
and the diffusion components respect ively. The relation between the surface 
potential <i>s, the gate-to-bulk voltage \ 'gb, and the source-to-bulk voltage Vsb 
can be given by ( 8 ] 


Vgb = Vfb + <t>s + lef/A 


(2.15) 


with [8] 

/l = yf#; s*(l - + V t h(e~' t,s / Vt>t + e^ s ~ <t ‘ , ~ VsB ^ v,h — e~^ i+VsB ^ Vth — 1), 

(2.16) 


where *) e ff is the effective value of the body-effect coefficient, and is a function 
of the channel doping and the oxide thickness, and <pi is the surface potential 
at the onset of strong inversion given by 


A = 2Vfc ln(N A /m) y (2.17) 

where N,\ is the substrate doping and is the intrinsic carrier concentration. 
In this model. 7 ,// has been expressed as [ 8 ] 

2 2 qUiNcjf jrsiLi.xiS fS '*-**nsiWL*w™' WK ) /o r«N 

%// = —77TT'~ * * 
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where C SILL, CSlLB , CSlWL , and CS1WE are fitting parameters neces- 
sary to model 7 <-// over a wide range of geometries. This equation is solved 
using an iterative numerical approach similar to the Newton-Raphson algo- 
rithm. Once <f>s is determined, the depletion charge Q' bS and the normalized 
potential F,s are obtained by the following two equations 


Qts = 7 tfjC' 0X \J <j>s + V^e' 

- Vm, 

(2.19) 

F iS = Vgb - Vfb ~ <t>s - 


(2.20) 


These values are then substituted for Q r b and F* in order to get E av at the 
source [Eqn.(2.8)J. This result is then used in Eqns.(2.6) and (2.7) in order 
to calculate the value of the effective mobility p e //s- In order to evaluate the 
surface potential <j>o and the related quantities such as the inversion charge, 
etc., at the drain end of the ohmic region (i.e., at y = !«„), the effective 
drain voltage Vox at this location must be determined first. If terminal biases 
are such that the device operates in the traditional linear region, then Vox 
is essentially equal to Vo, ignoring series resistance effects. On the other 
hand, in the saturation region, Vox is pinned to the saturation voltage Vo 3a t, 
which represents the electrical boundary between the linear and the saturation 
regions. Vo 9a t is the value of Vo at which the mobile carrier velocity attains the 
thermal saturation velocity, equal to v satl at the pinch-off point of the channel. 

With a fixed gate bias V G , as the drain voltage is increased, the inversion 
charge density Q\ D at the drain decreases, and the potential across the channel 
adjusts in order to increase the carrier velocity at the drain end so that the cur- 
rent flow is maintained. There is a limit as to how far Q' iD can decrease, since 
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the maximum carrier drift velocity possible is equal to the thermal saturation 
velocity v aai . 

Once Vo is large enough, so that the maximum carrier drift velocity v 3at 
is reached, Q' iD gets pinned at its minimum value, and the surface potential 
at the drain gets pinned at its maximum value 4>Dmax- The drain voltage 
corresponding to this potential is equal to Vo sat . As Vo increases beyond 
Vd* «t» the region C, shown in Fig.2.1, is formed. 

The excess drain voltage AV D (= V D — Vo 3at ) is dropped across this 
velocity saturated region. The length of this region l aat o increases with AVb, 
resulting in a finite output resistance. Vp 3at is obtained by equating the carrier 
velocity v at the drain end to be equal to v, at . Thus [8] 


v = 


[d 


( 2 . 21 ) 


where Fw is the normalized inversion potential at the drain end, given by 


F iD = Q'w/CL- (2-22) 

The expression for the drain current is given by Eqn.(2.13). In order to obtain 
a closed form solution for V Dsat , the following approximations are made in 
Eqn.(2.13). The square root term in the denominator of Eqn.(2.13) is assumed 
to be equal to unity, and faffavg is replaced by an effective mobility at the 
source end, faffs- These approximations are valid since the potential drop 
across the ohmic region of the channel is usually not very high. At Vd — Vp 3a t, 
4>o = <t>Dmax and the factor L re d is exactly equal to unity since the length of 
the velocity saturated region is still equal to zero. The carrier velocity at the 


15 



drain end can he expressed as (8] 


with 


v - v sat = 


B 


F\OminC, 


OX 


(2.23) 


B = + + (2.24) 


where is the minimum normalized potential required at the drain end 

in order to sustain the current Jo ■ Fuj m i n is related to <f>oma * by the feUewittg 
expression [8] 


FiOmtn = V(iB ~ Vf& ~ <t>Dmax ~ 7e//\/Wmox + - Vth- (2.25) 

Solving Etjn.(2.23), one gets a quartic equation for <f>omax of the following form 

[ 8 ] 


a $Dmax + b( t>Drnax + C^omax + d(p Dmax + e = 0, (2.26) 

where the coefficients are given by the following equations [8] 


a = a\, (2.2 6.a) 

b = — (2«j6i + 7e// a l)y (2.26 .6) 

c = b 2 + 2ajCi 4- 2 ')l u a x C 2 + Vt/i7e// a i> (2.26.c) 

d — — (2£qcq + d" 21 th 7 tff a iC 2 )i (2.26. d) 
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ami 


<’ «<■? + ' 'thl 2 '! f c%, (2.26.e) 

== l*tf fS j (2Li' 40 {), (2.26./) 

^1 = U.v//«i + <^5«i+2\ t/,ai + Fisai, (2.26.<j) 

n “ ^(.v// + lc.v//faat + 21'^VG e //a l + Fis<j>sdi~ 2VthFis<ii, (2.26./i) 

C 2 = 1 + ^ 50 i + 2VthO.\, (2.26 .z) 


V r 6V// = - Vfb- (2.26. j) 

These equations arc solved numerically and the root lying between fts and 
4>s + Yjjh is chosen, which gives the value of <j>D max , the maximum possible 
value of the surface potential at the drain end. Obviously, the value of 4>Dmax 
must lie between the surface potential at the source end <f>s and <f>s + Vdb, 
where \))n is the drain-substrate voltage. 

hi order to find Vusat, the following equation is used [8] 

* (2.27) 

where Yd b sat = Yp sat - Vb- The assumption that <pDmax is much larger than 
Vo, may not be satisfied in weak inversion. However, in this case, Vosat can 
be obtained very simply by noting the following two points. First, for Vos — 
V nStal (whore V DS sat - Yd sat - Vs), the drain current is given by [8] 


Id — Y r C' 0I Fi£)V sat , 


( 2 . 28 ) 



anti also, the drain current at the subthreshold region can be expressed by [11] 


Id = -jrV'fiVthC' OI F iS {l - e - v Ds.at/v th )^ (2 .29) 

The ratio F t i>( F,s is given as e~ VoSn, / Vth [11]. Using these relations (i.o., 
Eqns.(2.28) and (2.29)), V DSsat is obtained as V DSsat = V th In [1 +( v satL)/(ix e /fV th )], 
from which Yu, ai can be obtained by adding to it. With the saturation volt- 
age Vt)*at being known, the effective drain voltage Vbi can be found from the 
following equation [8j 

’ Dl = (i + {[Vo-v s )/v D „ t y-yi^ + Vs ’ (2 ' 30) 

where k is an empirical fitting parameter [8]. After Vox is determined, all the 
other drain side quantities can be found in a manner analogous to that applied 
for the source side. The drain end surface potential can be obtained by solving 
the following equation [8] 


VgB = V FB + 4>D + le/fA 


(2.31) 


where 


A = y/<j>£>(l — e~' i 'I Vth ) + V t h{e~*i>l v th 4- e (<f>D-<t>i-VDiB)/Vth _ e -Wi+vb lB )/vih _ i) } 

(2.32) 

where Vdib is given by Vdib = Vdi ~ Vi ?• The drain end depletion charge per 
unit area Q‘ hl) can be computed from the following expression [8] 



and the normalized potential at the drain end Fu > can be given by [8J 

F,» = Van ~ Vpb - <f>o - Q'bi>/C 0 x ■ (2.34) 

The effective mobility at the drain end Ht/fD can then be calculated by 
substituting the expressions for Q' bt> and Fio (Eqns.(2.33) and (2.34) respec- 
tively] for Q' b and F t respectively in Eqn.(2.8), and using Eqns.(2.6) and (2.7). 
The effective average mobility fi e /javg can then be calculated by substituting 
Hrffs and (i e no in Eqn.(2.14). 

The length of the saturated portion of the channel l, a tD, as reported in 
Fig.2.1, is obtained from a pseudo-two-dimensional analysis of this region as 
shown by El-Mansy and Boothroyd (12]. The solution of the two-dimensional 
Poisson’s equation for the surface potential <j> y in the velocity saturated region 
gives (12] 


<j> y = <l>D + FiatDceff sinh( - — ™ ). (2.35) 

where £>«.,/ / is the characteristic length which approximately represents the 
depth over which the mobile charge is spread in this region, and E aa t represents 
the lateral electric field at the boundary between the ohmic and the saturated 
regions of the channel, which can be expressed as E aat = rw/po, with /r 0 being 
the low- field mobility. At the drain end, i.e., at y = L, the surface potential 
can be given by 


<f>y = <t>D + {^D ~ t in)- 


(2.3(3) 
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Substituting Kqn.(2.36) in Eqn.(2.35), we get expression for L red , given by 


Lrrd = 1 — 


£cefj 

L 


sinh x ( 


YdjzVdi 

EsatDceff 


(2.37) 


Replacing the inverse hyperbolic sine function by its logarithmic equivalent and 
by using the approximation tanh(x) « x for x « 1, Eqn.(2.37) gets modified 
to (8] 

Lrrd = 1 - tanh[ ln(a + \/l + a 2 ) Dct J^ L ], (2.38) 


where a = ( V& — Vo^/iEsatDcefj). Equation (2.38) is of the form 1 — tanh(x), 
which can be simplified as 


1 - tanh(x) = 1 


e x - e~ x _ 2 

e x + e -x ~ e 2x + i * 


Therefore, Eqn.(2.38) can be reduced to 


(2.39) 


2 

red (a + \/l+ a i ) 2D “t// L + 1 

The parameter D ce j/ is known to be a function of the gate oxide thickness, 
doping density, drain junction depth, and applied bias [8]. Larger the gate 
field, the more strongly the carriers are accumulated near the surface resulting 
in a reduction in A*//- In this model, £><*// has been expressed empirically 
as (8] 



LCL 

1 + LCl(F iS - F w )Q' bS /C' 0X ' 


(2.41) 
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when* # 4I is the perrnitivitty for Si and LCL and LC 1 are empirical model 
parameters having the units of cm l/ 2 and V~ 2 respectively. 

In small geometry devices, it is important to model the Drain-Induced Bar- 
rier Lowering (DIBL) effect, which causes the surface potential at the source 
to he altered due to its proximity to the drain-substrate junction. Similarly, 
the surface potential at the drain is also affected by the state of the source- 
substrate junction. It causes the threshold voltage of a MOSFET to decrease 
as the reverse bias on the drain is increased. This results in an increase in 
the drain saturation current and a reduction in the output resistance as the 
drain bias is increased. In existing threshold voltage models, the DIBL effect 
is accounted for by making the threshold voltage V T a function of Vos [8], 
related by the following equation 


V T = Vtho - oV DS , (2.42) 

where Vtho it* the threshold voltage in the absence of the DIBL effect, and a 
is a parameter that controls the extent to which Vt is modulated by Vos- 

The present model accounts for the DIBL effect in the following manner. 
The surface potentials at the source and the drain ends are computed in the 
same manner as given earlier with the exception that Vqb is replaced by Vgb + 
A'anKob and by Vgb + KdibiVsB respectively, where K dib i is a model parameter. 
The DIBL effect cannot be modeled accurately using a single value for K dM 
that is adequate for all values of channel length L, since this effect increases 
exponentially with decreasing channel length [13]. In the present work, l\ dM 
has been modeled as [8] 

I< m = DIBLL x L dible , (2.43) 
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where I) I 11 L L arid D1BLE arc model parameters, and depend on the partic- 
ular device geometry. 

The simulation results and their comparison with the experimentally ob- 
tained results are presented in the next section. 


2.3 Analysis and Results 


The simulated model results are compared with the measured data available in 
literature for long channel ( L = 25 /im) [8] and for short channel (L = 0.55 /on) 
MOSFETs [8]. 

The simulated drain current (Id) versus the drain-to-source voltage (Vos) 
characteristics obtained from the model for different values of the gate-to- 
souree voltage {Vos) are shown in Figs.2.2 and 2.3 for the body-to-source 
voltage Vus = 0 V for two MOSFETs having channel lengths of 0.55 j um and 
25 fun respectively. The experimental data reported by Joardar et al. [8] are 
also shown for comparison. The results show a good match between the two. 
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Figure 2.2: //>4W characteristics obtained from our simulation for Vns = 0V 
and different values of Vos for a MOSFET having channel length of 0.55 f.nn 
and channel width of 25 ftm. The experimentally obtained results [8] are also 
shown for the sake of comparison. 
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Figure 2.3: I D -V DS characteristics obtained from our simulation for V BS = 0V 
and different values of V os for a MOSFET having channel length of 25 /mi 

an<l channel width of 25 fim. The experimentally obtained results [8] are also 
shown. 


ihe drain current (Id) versus the gate-to-source voltage (Vcs) character- 
istics for different values of the body-to-source voltage (V BS ) are shown in 
Figs.2.4 and 2.5 for the drain-to-source voltage \' DS = 0.5 V for two MOS- 
FETs having channel lengths of 0.55 fan and 25 fim respectively. Also shown 
in these figure are the experimental results obtained by Joardar et al. [8]. Our 
simulated model results match well with the experimental results. 


24 





Figure 2.4: IdA'gs characteristics obtained from our simulation for V DS = 
0.5 V and different values of V BS for a MOSFET having channel length of 
0.55 fim and channel width of 25 /mi. The experimentally obtained results [8] 
are also shown for the sake of comparison. 
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Figure 2.5: 7 0 -l bj characteristics obtained r 

0-5 and different vaincs of V f , “ 0 “ r Si “ Ulati0n f ° r = 

25 - and channe, W idth of 25 “ ^ *— - 

are also si,,™,, f ()r sake „ f Panmentaliy obtained results [8] 


‘eristics are sh^T ' ^ (%,) * arac ' 

- - - 

channel lengths of 0.55 m and 25 " ^ ^ 

obtained from Joardar at nf. (8J are also 

znr 
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Figure 2.6: g m -\ as characteristics obtained from our simulation for V B s = 0V 
and V DS = 0.5 V, for a MOSFET having channel length of 0.55 pim and channel 
width of 25 fim. The experimentally obtained results [8] are also shown for 
the sake of comparison. 
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Figure 1 2.7: g m -Vos characteristics obtained from our simulation for Vbs = 0 V 
and Vps — 0.5 V, for a MOSFET having channel length of 25 pm and channel 
width of 25 pm. The experimentally obtained results [8] are also shown for 
the sake of comparison. 


'Fhe output conductance (gts) versus the drain-to-source voltage {Vds) 
characteristics for different values of the gate-to-source voltage (Vos) are shown 
in Figs.2.8 and 2.9 for the body-to-souree voltage Vbs = 0 V for two MOSFETs 
having channel lengths of 0.55 pm and 25 pm respectively. The experimen- 
tal results reported by Joardar et al [8] are also shown in the figure. Our 
model results match well with the experimental results, and no discontinuity 
is observed in the characteristics. 
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Figure 2.8: i '««,• characteristics obtained from our simulation for different 

values of \'cs for Yds = 0 V for a MOSFET having channel length of 0.55 /rm 
and channel width of 25 fim. The experimentally obtained results [8] are are 
also shown for the sake of comparison. 


% 
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V DS (Volts) 

Figure* 2.9: gdr^us characteristics obtained from our simulation for different 
values of \ as for \ ’ns ~ 0 V for a MOSFET having channel length of 25 /im 
and channel width of 25 ftm. The experimentally obtained results [8] are also 
shown and compared with our model results. 


The MOSFET model developed by Joardar et al. was simulated and the 
measured as well as the simulated data are shown for a long channel and a 
short channel MOSFET. There is no discontinuity observed in the simulated 
drain current and conductance characteristics with respect to the terminal 
voltages. Also, a good match is obtained on comparison of model simulation 
results with the measured ones reported by Joardar et al. [8]. 


% 
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Chapter 3 


The Substrate Current Model 


3.1 Introduction 

With ever decreasing device dimensions in order to achieve faster and denser 
MOSFET lCs, the importance of proper modeling of the substrate leakage cur- 
rent has grown manifold. When a MOSFET is biased in the saturation region, 
a measurable substrate current I nb always flows. This current is the result 
of the electron-hole pair generation in the surface depletion region near the 
drain due to impact ionization. It is extremely sensitive to device dimensions, 
channel doping, terminal voltages, and temperature. Thus, there exists a need 
to accurately model the substrate current for proper circuit simulation. 

Several analytical substrate current models have been developed and re- 
ported in literature [14-16]. The substrate current model proposed by Aroxa 
and Sharnm [6] models the carrier ionization length l d near the drain in a semi- 
empirical manner. Since the substrate current /j U & depends exponentially on 
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this ionization length /*, a small error in its calculation can lead to a large 
error in the computed value of As Id is modeled semi-empirically, there 
is no physical significance for the model parameters, ^i> and J 2 appearing in 
101 - 


lu this work, a new physics based model for the ionization length Id is 
proposed, along the lines adopted by Wong and Poon [7]. This model for 
Id is then used in order to obtain an improved analytical expression for the 
substrate current /, u j. 


3.2 Model Formulation 

V W 

Cwwidort he MOSFET structure in saturation, z0 shown in Fig.3.1^ 

t ox represents the gate-oxide thickness, Id is the length of the velocity 
saturation region ( the impact ionization region), r.j is the drain junction depth, 
L is the effective channel length, and Vo, Vc, and V s represent the drain, gate, 
, ami source voltages respectively.^* 


-X fao su bs trat a curroa t I 3u b in an n-channel MOSFET can be given by [6] 




fC 


rL 

/ sub ^ Ip I dy , 

JL~Ia 


(3.1) 


where |f D is the drain currentfaricl a„ is the impact ionization rate of thelhigh. 
energy ) carriers infcm'Jlin the velocity saturated part of the channel (i.e., lj).^ 
^he^proximat^pression for thg^tfact ionization rate p^can begpieB toy-, 


** * 

r*' ? 
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o„ = .4, (T s ' /e , 


(3.2) 


where .4, and /i, art* ionization constants in cm -1 and V/cxn respectively. 
Therefore, the current expression [Eqn.(3.1)] can be rewritten as 

hub « In A, [ L e~ B ' /E(v) dy. (3.3) 

J L—ld 

For a MOS tiansistor operating in the saturation mode, the region near the 
chain is often called the velocity saturation region. The electric field in this 
region increases significantly and causes impact ionization due to high energy 
carriers traveling through this region. In order to obtain an expression for the- 
electric field distribution E(y) in the channel, consider Fig.3.1. 



Figure 3.1: Figure showing the coordinate system, velocity saturation region, 
ami boundary conditions, for evaluating the channel electric field E(y). 
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hollowing the analyses reported by Sonoda et al. [17] and Wong and Poon 
[i], we apply Gauss' law to the shaded region in the velocity saturated part of 
the channel, as shown in Fig.3.1, in order to obtain the following expression 

17 ] 


-E $at x } + E(y)xj + ^ f E ox (y) dy = ^ x j (y - L + l d ), (3.4) 

t$i JL-l d €$i 

where E OI (y) is the electric field normal to the Si/SiOz interface given by 
E„r(y) = V'r is the threshold voltage, V(y) is the potential along 

the channel, E sat is the minimum electric field for velocity saturation, and Na 
is the substrate doping density. 

In Eqn.(3.4), the electric field flux flowing from the velocity saturation 
region to the substrate is assumed to be negligible. The electric field in the 
channel E(y) is assumed to be independent of the depth x [17], and it is 
assumed to vary linearly from zero (at y — 0) to E aa t (at y = L — l d ). Thus, 
dE(y)/dy at y — L — l d can be approximated by [17] 


dE(y) , _ E sat 

dy ly ~ L ~' d L-l d - 


(3.5) 


Using the boundary condition given by Eqn.(3.5), and following the approach 
similar to that adopted by Sonoda et al. [17], one gets the following equation 
for the channel potential V (y) 


\ 2 p 

V(y) = (cosh (D) - 1) 4- XE sat sinh (D) 


where V Dsat is the drain saturation voltage at y = L - l d , D is equal to 
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(.V - !■> + frf)/A, anti A, the characteristic ionization length, is defined as [7] 


A = J'f-xj (3.7) 

y c ox 

Differentiating Eqn.(3.6) with respect to y, one obtains the expression for the 
electric field E(y ), given by 

£(y) = -z — -r- sinh (D) + E, at cosh (D). (3.8) 

L — l<i 


From Ecjns.(3.6) and (3.8) and using the identity cosh 2 (D) — sinh 2 (D) = 1, one 
gets the following relation for E(y) 


£(y) = 


1 


/ YSmI Yjoibi 4_ 

1 A L-l d 


) 2 + El 


X 2 El 


sat 


sat 


(L - Uf 


(3.9) 


Differentiating Eqn.(3.9) with respect to y, one gets 


iy = . xdE (3.10) 

Jew 

Changing the variable y in the expression for (Eqn.(3.3)] to E one obtains 


I sub 


MX f 

JE 


e -BiJE(y) 


Jm 2 - Mu+v 


dE, 


(3.11) 


where b is equal to A E tat /(L - l d ), E sat is the electric field at y ~L-l d , and 
E m is the maximum electric field at y — L. 
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iho integral appearing in Eqn.(3.1 1) does not have a closed form solution. 
However, based on the assumption that E m is much greater than E sat , the 
expression for J 3ub can be approximated as [6] 


/ 


.1**6 



(3.12) 


In h(jn.(3.12), the term E] at and b 2 can be neglected as compared to E^, and 
the simpliiied expression for I sub can be given by 


Hi 


(3.13) 


In this work, the maximum electric field E m has been modeled empirically 

as 

r _ V DS ~ V^Daat 

m — m — • (3 - 14) 

where t/ and 0 are dimensionless fitting parameters. Thus, substituting the 
expression for E m (Eqn.(3.14)] into the expression for J 3ub [Eqn.(3.13)], a new 
analytical expression for the substrate current I 3ub has been obtained in this 
work, which is given by 


= Ijo^Vds - n Vo,.') E e -W»./C‘W-,W). (3. 15 ) 

The ionization length Id can be determined as a function of the channel length, 
A, and Vp - V Daat (7). In order to determine this relation, consider Eqn.(3.6), 
which gives the potential as a function of the position in the channel. 


3G 



Al y /., \ (»/) is equal to the drain voltage l )). Substituting this in 
Kt|ii.(3.(>)> we get the following expression 

= r^ (cosh( i) _ [) +sinh ( i>- (»•«) 

Substituting X/(L - Id) as a, (Vp — Vp 3at ) / (XE 3at ) as b, and e ld ^ x as c in 
K<|ti.(3.1G), one obtains the following relation 

2/> = «(r+ - -2) + (r- -). (3.17) 

c c 


It is to be noted that all the three terms a, 6, and c depend on Id- However, 
since the term c is exponentially dependent on Id, hence, its variation with U 
would be much more drastic than the variations in the terms a and b with the 
variation in /<*. Therefore, Eqn.(3.17) is expressed as a quadratic equation in 
c and the solution can be given by 



(3.18) 
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whi’ii* A is a ilinu'iisionlrss model parameter which depends on Vp — Vp, at . In 
this wotk. A" has been modeled as 


A — A|(\y — Vp sa t) 2 + K?(Vp — Vp aat ) + K 3 , (3-21) 

where k\ {V' 2 ), h'x (l' -1 ), and A3 (dimensionless) are fitting parameters. 

In the evaluation of the substrate current I snb , Vp, Q( is an important pa- 
rameter which must be determined. It can be modeled as [6] 

Vom = + - l), (3.22) 


where o is a fitting parameter. 

The improved model developed for I ittb in this section was simulated and 
the results obtained were compared with the experimental data given in [6]. 
The results are presented in the next section. 


3.3 Analysis and results 

For the calculation of the values of the substrate current I 3u b as a function of 
Vos, the parameters required to be extracted are Ap Bp 77 , A, /*, P, and a 
(Eqn.(3.15)]. The device parameters used for the simulation are Xj = 0.4 firn 
and t ox = 150 A 0 , which gives the value of A to be equal to 1.37 fxm [from 
Eqn.{3.7)}. 

In this work, the simulated results obtained from the model developed for 
the ionization length l d (Eqn.(3.19)j are plotted as a function of V D - V D , at 
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m F ig 3 2. 'Hm* symbols in Fig.3.2 are the data generated from MIN1MOS, 
ubi aim'll fiom Wong ami Poon [7], for the channel length L = 0.5 pm. From 
a comparison of the two, the values of the fitting parameters Ki, K^, and K 3 
are obtained as -0.03 V -2 , 0.27 V -1 , and 1.2 respectively. 



Figure 3.2: ionization length as a function of V D - V D3at) for L = 0.5 pm. 
Symbols indicate MINIMOS generated data [7]. 

The values of the ionization length Id obtained from Eqn.(3.19) are also 
plotted as a function of the channel length L in Fig.3.3, for Vp — Vg > aot = 2.0 V 
and 2.5 V. The data generated from MINIMOS [7] are also shown in the figure. 
The results show a good match between the two. 

The values of the fitting parameters K u and K z extracted above are 
used to determine the ionization length Id- The values of the ionization con- 
stants A t and B, are kept equal to 2.54 x 10 6 cm -1 and 1.92 x 10 6 V/cm re- 
spectively, and those of the saturation electric field E 3a t and o are kept equal 
to 1.053 x 10 1 V/cm and 1.19 respectively, as given by Arora and Sharma [6]. 
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Figure 3.3: Ionization length I 4 as a function of the channel length L, for 
l/j - I that = 2.0 V and 2.5 V. Symbols indicate MINIMOS generated data 


hi order to extract the values of the fitting parameters /? and 77 used in 
our empirical expression for the maximum electric field E m [Eqn.(3.14)], the 
exact equation for I tU b [Eqn.(3.11)] is integrated numerically and the results 
are plotted in Fig.3.4 as a function of Vos- Then, the values obtained from 
the approximate equation for I su 5 (Eqn.(3.15)] are mapped onto it, and are 
also shown in Fig.3.4. From this mapping, the values of fi and 7 ; are obtained. 
Table 3.1 shows the values of the extracted parameters needed for determining 
the substrate current /,„ j. 
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Figure 3.4: The values of /,„& obtained from the approximate relation 
(Eqn.(3.15)j, and the numerically integrated values obtained from the exact 
expression [Eqn.(3,ll)] versus the gate-to-source voltage V C s, for three differ- 
ent channel lengths (0.5 pm, 0.6 pm, and 0.7 pm). 


Table 3.1: The values of the parameters needed for determining the substrate 
current 


A, 

io 6 (o/r 1 ) 

Bi 

10 6 (V/cm) 

E ta t 

10 4 (V/cm) 

V 

a 

0 

2.54 

1.92 

1.053 

1.32 

1.19 

0.244 


The simulated results obtained from the analytical model equation for I 3ub 
[Eqn.{3.15)J, using the values of the model parameters extracted earlier, is 
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compared with thr experimental data given by Arora and Sharma [6]. Figure 
3 5 slmvvs tin* substrate current I mb obtained from our model as a function of 
the gate- to- source voltage \ for Vps - 4.6 V and Vus = 0 V. The cxperi- 
ineiital data reported by Arora and Sharma [6] is also shown for comparison. 
Our model results match well with the reported experimental results. 



Figure 3.5: The simulated values of the substrate current hub obtained from 
our model as a function of the gate-to-source voltage V GS , for V DS = 4.6 V 
and Vus * 0 V. The experimental data reported by Arora and Sharma [6] 
are also shown for comparison. 

The nature of the substrate current is such that it initially increases with 
increasing V as and then it decreases as V GS increases beyond a certain value, 
thus showing a peak at a particular value of \ os- This can be physically 

explained as follows. 

The substrate current hub is a function of Ids, V DS -V DsaU and l d [Eqn.(3.15)J. 
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Initially, f* *r a small \'<;s, Ids increases rapidly with increasing Vast thus Liub 
hh-U’um's atmtdtngly. However, as \\;s increases, Vos ~ V o sa t and Id both 
deneiiw* Heme, a point is reached when the effect of increasing I S ub due 
to increasing \\;s is nullified by the decrease in I su b due to the decreasing 
\ ), s - \ ) Ha! and la- Then onwards I sub decreases as V C s is increased. 

A simple and accurate approximation of the ionization length Id, based on 
the calculation of the electric field distribution near the drain region has been 
developed in this work. Our model results are verified with the MINIMOS 
simulation results reported by Wong and Poon {7], and a good correlation 
between the two is obtained. 

This new model for Id is then used to calculate the substrate current. The 
results obtained from the improved substrate current model are then compared 
with the experimental data reported by Arora and Sharma [6], and a good 
match is obtained between the two. 



Chapter 4 


Transconductance Overshoot 
Model for Submicron MOSFETs 


4.1 Introduction 

Different analytical submicron MOSFET models have been developed in 
order to incorporate the various short channel effects. Most of these models 
accurately reproduce electrical characteristics of short channel devices down 
to half micron channel lengths. However, they encounter some difficulties in 
accurately reproducing electrical properties for devices with channel lengths 
less than 0.1 fim. 

One such effect that becomes predominant when device dimensions shrink 
is the electron velocity overshoot. It is one of the most important effects from 
the practical point of view, as it is directly related to the increase in the current 
drive and the transconductance of the device, which have been experimentally 
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ultM'i vcd hi shml channel MOSFIvTs. Some authors (1-4) have shown that the 
experimentally measured values for the transeonductance for submicron MOS- 
1'K l s (fur channel lengths less than 0.15 /im) are higher than the theoretically 
predicted maximum values that can be reached. 

Figures 4.1 ami 4,2 show the drift velocity of the electrons as a function of 
time and distance from the source respectively in silicon for different values of 
the electric Held jl8j. From these two figures, it can be seen that for electric 
(Mils greater than 10 kV/cm, there is a possibility of velocity overshoot for 
electrons if the channel lengths are short enough. 



^Figure 4.1: Variation of the drift velocity of electrons with time, for different 
values of the electric field E [18]. 

•IV principal aim of this work it to develop a model for the trauscomh.c- 
tance overshoot observed in ultra-short channel MOSFETs. In this work, the 

carriers consideration are electrons. The MOSFET transeonductance 

( s ,„) is the most important figure of merit when dealing with large signal 
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Distance from source (Jim) 


Figure 4,2: Variation of the drift velocity of electrons as a function of distance 
from the source for different values of the electric field E [18]. 

switching performance of logic devices, as the time constant for MOSFETs to 
charge a load capacitance C is proportional to C/g m . This parameter is also 
of immense importance from the analog circuit point of view, since the gain of 
tlu* MOSFET amplifiers is directly proportional to g m . 

Many models have been developed and presented in literature, which in- 
corporate the carrier velocity overshoot effects [5,19,20]. The drift velocity v 
of electrons in an inhomogeneous electric field has been expressed as [5,21] 



where A is the parameter proposed by Thornber [21], which depends on the 
longitudinal electric field, tty is the drift velocity in a homogeneous field, and E 
is the electric field in the channel. The corresponding drain current expression 
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is giu*ii In [V 


/„ - WfoWaS' Vits) t II' 

/(I + taitisT - + ^ ds), (4.2) 

Ivlin.* A. IS l III- vi'lurity oveishoot parameter (rm 3 /V-sec) as given by lloklan 


rf tt/ ^ is t,H * km ' fit,,d c!ectro » mobility, u, a< is the saturation velocity of 
elect runs, ir IS I lie i Imuiicl width, L is the channel length, and 


f : (V<;s,\'i>s) = f' DS Q'(V)d\ 


(4.3) 


when* </(l ) is the mobile channel charge per unit area. 

l*rom hqn.(4.2), the MOSFET transconductance g m can be given as [5] 


!hn ~ 


#fo . a.. 

L 0l' 6 . s l i+ + L J ' 


(4.4) 


This expression assumes that the low-field mobility /i 0 is not a function of 
the applied gatc-to-sourcc voltage Vqs- This is contrary to the conventional 
mobility model (11), which accounts for the dependence of mobility on Vqs, 
given by 


jjQ 

= r+eivcs'-ti-y (4 ' 5) 

where ft (! is the mobility as a function of the gnte-to-source voltage Vas, Vr 
is the threshold voltage, 6 is the mobility degradation coefficient [11], and /x 0 
is the low field mobility at \ };s - Vr- Tsividis [11] has given the empirical 
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I’\ pi «'*»**H Ml f* »r 0 .|.s 


0 = A (4.6) 

tox 

when* f Im* »*r <ig typically varies from 0.001 to 0.004 pm/V [11], and t 0 * 

is i in' oxide 1 i»*i kn***»s Alsu, the parameter A a has been shown to be functions 
uf the low- held mobility and the temperature [5], however, the dependence of 
this parameter on the bias voltage is not shown. 

The aim of this work is to develop an analytical mobility model in order to 
explain the phenomenon of higher than theoretically predicted transconduc- 
tance observed for ultra-short channel length MOSFETs. 


4.2 Model Formulation 

Assume that the source is located at x = 0 and the potential at that point 
is zero, The drain is assumed to be at x = L (where L is the effective channel 
length) and is at a potential V D s- The drain current Io can be calculated as 
a product of the mobile channel charge per unit area, the channel width, and 
the electron velocity, i.e., the electrons in transit model, given by 

Id « Q'{V)vW. < 47) 

The illifi velocity mode! used in this work is given by [5] 

ftE (4.8) 

v “ TT~K' 
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uhnr /: IS thi* electric field in the channel, E c us the critical electric field 
fur 'i lunty saturation, ami p is t he enhanced electron mobility, a model for 
ttlmh is |iinjn »sis I later in this chapter. This drift velocity model predicts 
an almost linear increase of the velocity v with the electric field E until the 
critical elect tie field • is reached, beyond which the velocity saturates to the 
sanitation velocity r,„|. This behavior of the carrier velocity v as a function of 
the elect tic field E is shown in Fig.4.3. For getting the characteristics of Fig.4.3, 
the values uf and r (#l> which are substituted in Eqn.(4.8), are assumed 
to be equal to 25 k\'/nn and IQ 7 an/s respectively. Also, an approximate 
expression fur h'c can be given as Ec =* Vioe/^Gi where no is the mobility 
given by Kqn.f 1.5), 



Figure 4.3: Variation of the drift velocity of electrons with the electric field. 


The values used are Ec = 25 kV/cm and v, at = 
!* as filC ami for high fields, t> « v m t- 


10 7 an/s. For low fields, 


In this work, we propose a new model for the mobility as follows. Hie 
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t infill »>f Hie proposed equation is semi-empirical. Physically, it is obvious that 
thr mobility would imrca.se with larger applied drain-to-source voltage V D s 
ami .smaller channel length L. However, at the same time, it is a function of 
the critical electric field Ec- Thus, the proposed model is given as 


it = n c {l + 


Leo Yds 
Ea> L? 


). 


(4.9) 


where /<<; is mobility as given in Eqn.(4.5), and L co and £«> are empirical model 
paraiiieters Substituting Ktjns.(4.8) and (4.9) in Eqn.(4.7), we get 


h = wwn 


Hg(1 + 


1 + 


kaiYjO&\ 
EcoLX i t? 

V*at 


(4.10) 


Rearranging trim*. we got 




L m Vds . 


(4.11) 


In the term [1 + (/icE/e.*)] in Equ.(4.10), the electric field B can be assumed 
to have a constant value, which is equal to Vds/L. This assumption is valid 
since as E changes, there is only a slight variation in the term [1 + 0*7*7 «»<)] 
as. in general, po E/v.. t « l. Substituting E = V BS /L and integrating the 
left-hand side of Eqn.(4.11) from 0 to L and the right-hand side from 0 to Vos, 
in a manlier similar to that adopted by Roldan et al. [5], we get 


/; w , + ^>dr= /;“ nw)Mi + w < 412 > 


or 


50 



(4.13) 


, = «; {Mj + IgOTWjM 

L u eatog 

T v„,t 

hquainm 1 1 13) gives t hr expression for the drain current, which was derived 
using the new model for the mobility proposed in this work. The transcon- 
liuetiitire ran lie obtained by taking the derivative of Eqn.(4.13) with respect 
to \ "[ ; \ , whieh ran be given by 




W. , l, 

t (1 + e„ 


ro 1 ns 


L 2 




(-1. 1-1) 


when* 


. (<}FI<n\ ;s ) it<;+F{dtt«/d\cs) fF/i G 7 05 /(t; iat L)](d/x G /dK G s) . . 

’’ * ^ _j_ noVos Y ’ ' ' 


i 4 . euiia 


whirlt run be expressed as 


_ (tlf/fll 0 , 

i + f-iilui ' l + ea!iu ; ’ 

‘ r - -* ' V.alL 


(4.16) 


v,»iL 


when* J is given by 


& 


^^g[&Ygs) 

{0F/dV G s)tiG 


(4.17) 


Differentiating the expression for Ha (Eqn.(4.5)] with respect to V^s, we get 


dfic zM. (4.18) 

d\’as ~ l 1 + 0( v gs - It)] 2 
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S 1 1 1 •*> * 1 1 « 1 1 n 1 1 *, I qit ( 1 lit) in l\q».{4.14), we got the final expression for the 

t I.IIIM < ■)!■ 1 1 1< I .IlM V AS 


fhn 


u 

7. 


(i * 


l;.Vfs A0F/i n'cs)HG t . , l , 

K ta L 1 1 + + 1 + ugXdjl'' 


(4.19) 


V„tl 


V,at L 


Equation ( l ID) is th«’ oxprossion for the MGSFET transconductance, which in- 
corporates I ho improved mobility model proposed in this work. The transcon- 
tluc iant o mi’tshooi f.n tor I\ can be given as 


K = 1 + 


Leo I ns 

ico L* * 


(4.20) 


whirl) is obviously gt outer than unity, and, in turn, predicts a higher than 
expiM toil trnuseonduetnnoe. 


4.3 Results and Discussion 


The simulated values of the transconductance obtained from Eqn.(4.19) are 
compared with the experimental values for a MOSFET reported by Sai-Halasz 
i t id, |l) as a function of the channel length. This comparison is shown in 
Fig,.-*.-*, which shows the transconductance per channel width (g m /W) as a 
function of the channel length L for two different temperatures (77 K and 300 


K). In the figure, the experimental values [1] are denoted by symbols. From 
Fig. 4. 4, it can be seen that as the channel length decreases, there is a sharp 
increase in the transeondiictancc. Our analytical model results match quite 
well with the experimental values reported by Sai-Halasz et al. (1]. The value 
of the parameter W/i)i;« is kept the same as that given by Roldan et al. [5] 

hik! is given in Table 4.1. 
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Figure I I The trausrouduciance per channel width reported by Sai-Halasz 
ft al jl] as a fum lion of the channel length for two different temperatures 
(77 K ami iilKI h) arc shown by symbols. The results simulated from our 
model lliucj.) arc also shown for comparison. The parameter values used in 
the smmlntioii are \ - Vf - 0.6 l', Vps - 0.8 V, and t ox — 4.5 nm. The 

other model parameter values arc shown in Table 4.1. 

In the enlr ulatum fur the transconductance g m , the value for the low-field 
mobility ft a is assumed to be 390 rr« a /V- sec and 720 cm 2 /V- sec for T = 300 K 
and 77 K u-spertively |5) The value of jig [appearing in Kqn.(4.6)J is chosen to 
be equal to 0 002 pm/t and the value of the parameter 9 can be obtained from 
K<|ti ( f (>) for a given value of f M . The values of and 6 are then substituted 
in Kqn (4.5) in order to obtain the value for pc. The critical electric field E c 
is then computed from the expression Ec — v sa t/HG- 

‘f'he model parameters used in the expression for the transconductance 
[Kipi ( 1, 1*1) j a»e l 0, and (i. In this work, tlie values of the parameter 
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E,„ is taken to bo equal to that of the critical electric field Ec - This can be 
justified as follows. The increase in the transconductance is observed only for 
high electric fields, i.e., for E > Ec- In our model, the drift velocity satu- 
rates for E greater than Ec- Therefore, in order to account for an increase 
in the transconductance, the mobility must increase [Eqn.(4.19)]. This is ac- 
complished by keeping the parameter E co to be equal to Ec- The values for 
the other two fitting parameters L co and p are extracted from a match of 
the experimental and the simulated transconductance values. The extracted 
parameter values for this case are shown in Table 4.1. 

Table 4.1: Extracted model parameter values obtained from a comparison 
of the results simulated from our model with experimental transconductance 
values reported by Sai-Halasz et al. [1]. 


Temp. 

(A') 

ft o 

(rm 2 /V-. sec) 

f‘a 

(an 2 /V- sec) 

Eco 

(kV/cm) 

■f'ro 

(10 -6 cm) 

P 

V 3 at 

(10 7 cm/s) 

dF/dVcs 
(10~ 7 C/cm 2 ) 

300 

390 

307.9 

25.98 

0.935 

0.48 

0.8 

5.2 

77 

720 

568.4 

17.59 

0.791 

0.48 

1.0 

5.5 


The results obtained from our model are also compared with the ones re- 
ported by Pinto et al. [22], The intrinsic transconductance versus channel 
length data given by Pinto et al. [22] were calculated by using accurate sim- 
ulators which included the velocity overshoot effect. The simulator PADRE, 
used by Pinto ct. al. [22], made use of the energy transport equations [23-25], 
defined in terms of the average carrier drift velocities and temperature. 

Figure 4.5 shows the transconductance per channel width g m /\V obtained 
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from our model, as a function of the channel length L, for different values for 
the low-field mobility (125, 250, and 500 an 2 /V -sec) at T = 300 K. The results 
reported by Pinto et al. [22] are also shown in the same figure for the sake 
of comparison. The parameter values are extracted after matching our model 
results with the simulated transconductance values reported by Pinto et al. 
[22], and are shown in Table 4.2. 



Channel Length (pm) 

Figure 4.5: The simulated intrinsic transconductance per channel width re- 
ported by Pinto et al. [22] as a function of the channel length for three different 
values for the low-field mobility (125, 250, and 500 an 2 /V -sec) at T — 300 K, 
shown by symbols. The results simulated from our model (lines) are also shown 
for comparison. The parameter values used in the simulation are V G s = Vds = 
1.5 V, v sat = 10 7 cm/s, dF(V G s,V D s)/dV G s = 10 -6 C/cm 2 ,t ox = 4 nm, and 
\' T « 0.45 V. The other parameter values are shown in Table 4.2. 
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Tabic 4.2: Extracted parameter values after comparison of our results with the 
simulated transconductance values reported by Pinto et al. [22]. 


Ho (cm 2 /V- sec) 

Hg (cm 2 /V-sec) 

E co ( kV/cm ) 

L co (10 -6 cm) 

P 

125 

81.96 

122 

2.44 

0.48 

250 

163.9 

61 

2.26 

0.48 

500 

327.85 

30.5 

1.34 

0.48 


The transconductance overshoot factor F v computed from Eqn.(4.20) is 
plotted as a function of the channel length for two different temperatures (77 
K and 300 K) in Fig.4.C along with the experimentally reported values by Sai- 
Halasz et al. [1]. It can be seen from the figure that for large channel lengths, 
the value of the overshoot factor F v is close to unity, indicating that there is 
no increase in the transconductance. However, for channel lengths less than 
0.15 /mi, the overshoot factor becomes quite significant. Our model parameter 
values are extracted after matching Eqn-(4.20) with the experimental values 
reported by Sai-Halasz et al. [1], and are presented in Table 4.3. 


Table 4.3: Extracted model parameter values after comparison of the overshoot 
factor simulated from our model with the reported values by Sai-Halasz et al. 

(I)- 


Temp. (K) 

E C o {kV/cm) 

Leo (10 -6 cm) 

77 

0.791 

17.59 

300 

0.935 

25.98 


5C 























Channel Length (}im) 

Figure 4.C: The room temperature (300 K) and low temperature (77 K) over- 
shoot factor F v as a function of the channel length, as reported by Sai-Halasz 
et al. [1], shown by symbols. The results simulated from our model (lines) are 
also shown for comparison. The parameter values used in the simulation are 
Vos - Vr — 0.6 V and V DS — 0.8 V. The other parameter values are given in 
Table 4.3. 


The results obtained from our analytical transconductance model match 
well with the experimental data reported by Sai-Halasz et al. [1], as well as 
with the simulated transconductance behavior reported by Pinto et al. [22]. 
It has been observed in this work that the values of the model parameter L c0 
change with the low-field mobility //o, and, hence, also with /«<;. Table 4.4 
shows the variation of this parameter with /jq- 
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Figure 4.7: Extracted values of the parameter L co (shown by symbols) as a 
function of the mobility no along with a straight line best fit. 

model the variation of the mobility itself with the channel electric field. This 
mobility is also a function of the drain-to-source voltage, oxide thickness, gate- 
to-source voltage, and channel length. With large channel electric field and 
small channel length, the mobility is bound to increase, and the increase in 
the transconductance is then accounted for by the correspondingly enhanced 
electron mobility. 

Equation (4.20), which gives the expression for the overshoot component of 
the transconductance, is seen to be inversely proportional to the channel length 
and directly proportional to the electric field in the channel, which vindicates 
our earlier statement. Thus, E co and L co can be thought of as representing 
something similar to the critical values of the electric field and the channel 
length respectively. For electric fields comparable to or larger than E co , and 
for channel lengths comparable to or smaller than L co , the overshoot factor 
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cannot be neglected. 


The parameter L co is observed to be a function of the electron mobility 
fjta- In the analytical expression for L^, A and B have been used as fitting 
parameters due to the absence of any experimental data. However, the straight 
line fit is not accurate, as is obvious from Fig.4.7, and, thus, the extracted 
values of the parameters A and B may have some errors in them. 



Chapter 5 


Summary and Conclusions 


The growth in VLSI has given rise to the need for developing accurate 
MOSFET models, which are capable of explaining device behavior down to 
deep submicron channel lengths adequately. Therefore, accurate MOSFET 
modeling has proved to be an essential step towards design and analysis of 
electronic circuits. 

With device dimensions being scaled down, there is a corresponding in- 
crease in the drain-to-substrate leakage current 7 3u6 as a function of the drain- 
to-source voltage, resulting from the impact ionization of high energy carriers 
in the surface depletion region near the drain. This current is the result of the 
electron-hole pair generation in this region near the drain, and is extremely 
sensitive to device dimensions, channel doping, terminal voltages, and tem- 
perature. An immediate consequence of this substrate leakage current is the 
increase in the drain-to-source output conductance g da . This parameter is very 
important from the analog circuit point of view, as it directly affects the am- 
plifier gain, which uses these devices. Hence, an accurate model for 7 su f, is 
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necessary. 


Another effect that comes into prominence with decreasing channel lengths 
is an abrupt increase in the device transconductance g m for ultra-short channel 
length MOSFETs. Not only is this parameter of immense importance from 
the digital circuit point of view, but more so from the analog circuit point of 
view. In digital circuits, the time constant required for a MOSFET to charge 
a load capacitance C is proportional to CJg m - Thus, g m is an important factor 
in determining the speed of operation of digital circuits. In analog circuits, 
since the gain of the MOSFET amplifiers is directly proportional to g m , any 
change in this parameter affects the amplifier gain. Hence, this behavior needs 
to be modeled adequately. 

In this work, these two effects have been modeled in a semi-empirical man- 
ner. Prior to this, an existing physics based MOSFET model reported by 
Joardar et al. [8] was rederived and the results obtained from our simulation 
were compared with the experimental data reported [8]. A brief summary of 
this work is given below. 


5.1 Summary of the work 


A major drawback in the existing MOSFET models is the presence of 
discontinuities in the derivatives of the drain current with respect to the gate- 
to-source voltage and the drain-to-source voltage. These discontinuities arise 
primarily due to the nature of the velocity-field expressions and the smoothing 
functions used in the transition area between the linear and the saturation 
regions in these models. The model developed by Joardar et al. [8] attempts 
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to overcome the problems mentioned above. This new model is normalized 
surface potential based and is quasi-static in nature, i.e., the variation of the 
terminal voltages is assumed to be sufficiently slow, so that the they can be 
assumed to be identical to the dc voltages at any given time. This model was 
simulated and the results obtained were compared with the experimental data 
reported by Joardar et al. [8]. 

A simple and accurate approximation of the ionization length lj is (level- 
oped in this work. This has been obtained by a calculation of the electric field 
distribution near the drain region. Our model results were verified with the 
MINIMOS simulation results reported by Wong and Poon [7], and a good cor- 
relation between the two was obtained. This new model for Id was then used 
to calculate the substrate current. The results obtained from the improved 
substrate current model developed in this work was then compared with the 
experimental data reported by Arora and Sharma [C], and a good match was 
obtained between the two. 

In this work, the abrupt increase in the transconductance g m with decreas- 
ing channel lengths has been accounted for by proposing a new semi-empirical 
model for the electron mobility. The mobility has been made a function of the 
(irain-to-source voltage, oxide thickness, gate-to-souroe voltage, and channel 
length. With large channel electric field and small channel length, the carrier 
mobility is bound to increase, and the increase in the transconductance is then 
accounted for by the correspondingly enhanced electron mobility. Our model 
results have been compared with the experimental data reported by Sai-Halasz 
at al. [1] and the results showed a sufficiently good match. 
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5.2 Scope for Improvement 


The semi-empirical substrata current model developed in this work has 
expressed the maximum electric field E m in the channel in a semi-empirical 
manner, with 7/ and 0 being the fitting parameters. These parameters can 
be expressed physically as a function of the oxide thickness, substrate doping, 
channel length, and drain-to-source voltage. Also, there are several approxi- 
mations made in evaluating the electric field E(y) as a function of the position 
along the channel [I\qn.(3.8)j. The Gaussian box in Fig.3.1 is assumed to ex- 
tend only upto the drain-junction depth Xj, whereas in reality, the depletion 
region extends well beyond Xj. The electric field flux flowing from the velocity 
saturated region to the substrate is assumed to be negligible, and the electric 
field within the channel E(y) is assumed to be independent of the channel 
depth x. Thus, a scope for improvement exists in the accurate determination 
of the electric field within the channel. 

In the mobility model developed in this work, the carriers taken into consid- 
eration are electrons. The same results can be extended for obtaining the new 
enhanced hole mobility, and, subsequently, the transconductance expression 
for a p-channel MOSFET. In the expression for the electron mobility devel- 
oped in this work, some physical significance has been given to the parameter 
E co , but the exact physical significance of the parameter L co is not yet clear. 
Also, this parameter has been approximately expressed as linear function of 
the electron mobility yc in a totally empirical manner. The future scope for 
improving the mobility model developed in this work lies in imparting physical 
significance to the parameter L c0 , as well as obtaining a better model for L co 
as a function of ya- 
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